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We study numerically the eflFects of the Rashba spin-orbit interaction on the model of electrons in 
n-doped semiconductors of Matsubara and Toyozawa (MT). We focus on the analysis of the density 
of states (DOS) and the inverse participation ratio (IPR) of the spin-orbit perturbed states in the 
MT set of energy eigenstates in order to characterize the eigenstates with respect to their extended 
or localized nature. The finite sizes that we are able to consider necessitate an enhancement of the 
spin-orbit coupling strength in order to obtain a meaningful perturbation. The IPR and DOS are 
then studied as a function of the enhancement parameter. 
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The metal-insulator transition (MIT) is one of the 
paradigms of Condensed Matter Physics P, 0| and new 
features constantly appear according to the physical 
properties under study, the specific system or the emerg- 
ing experimental techniques. The richness of the physics 
around the MIT stems from the fact that it is a quan- 
tum phase transition where disorder and Coulomb inter- 
actions coexist and compete in the determination of the 
ground state. In the case of the n-doped semiconductors, 
the MIT appears at doping densities where the Fermi 
level is in the impurity band j3|J. This observation al- 
lows a description taking into account only the electronic 
states built from the hydrogenic ground state of the dop- 
ing impurities. For densities slightly larger than the 
critical one (in the metallic side of the transition) non- 
interacting models, like the Matsubara- Toyozawa (MT) 
Q are applicable. Furthermore, the previous description 
in terms of impurity sites is often traded by the Ander- 
son model of a tight-binding lattice with on-site or hop- 
ping disorder. Large amounts of numerical work have 
been devoted to the Anderson model [Hj and the criti- 
cal exponents obtained fit reasonably well those of the 
experimental measurements 

The recently developed field of spintronics is contribut- 
ing to put the MIT again into the focus of the condensed 
matter community. A key concept for possible applica- 
tions of spintronics is the spin relaxation time, that is, the 
typical time in which the electron spin loses its initially 
prepared direction. Interestingly, the maximum spin re- 
laxation times in n-doped semiconductors have been ob- 
served for impurity densities close to that of the MIT 
0-[lll. This intri guing physics is not completely under- 
stood at present, and various mechanisms of spin relax- 
ation have been thought to be active at the MIT region 
12l4l5l|. At the level of models, the generalization of the 
Anderson model in order to include some spin-orbit cou- 
pling has been provided by Ando [16[. While this model 
is very useful to study the progressive breaking of the 



spin symmetry 17|, its connection with experimentally 
relevant systems requires the estimation of coupling pa- 
rameters which are not obtainable from first principles. 
This situation has lead us to reconsider the MT model of 
impurity sites randomly placed in order to incorporate in 
it the spin-orbit interaction. The various spin-orbit cou- 
plings (intrinsic and extrinsic) can be included and lead 
to effective Hamiltonians which depend on fundamental 
material constants, rather than on adjustable parame- 
ters. 

In this paper we first consider the MT model in order to 
characterize the regions of extended and localized states, 
analyzing the limitations of the model and the conditions 
of applicability. We then include one of the sources of 
spin-orbit coupling, the Rashba interaction [H] to study 
how the previous picture evolves under increasing val- 
ues of the coupling strength. This work is a necessary 
step towards the understanding of spin dynamics in the 
generalized models that will allow us to extract the spin 
relaxation times close to the MIT. 

We start by considering the MT Hamiltonian Q 
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where 



(crncr) represents the creation (annihilation) 



operator for the ground state of the impurity at site m' 
(m) with spin projection a in the z-direction. The spin 
degree of freedom is irrelevant for the MT model, but it 
will become crucial later. The hopping matrix element is 
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where 0p(r) = (j){\r - Rp|), with (/)(r) = l/{Tra^y/'^x 
exp(— r/a), and a is the effective Bohr radius. The 
Coulombic potential produced by the impurity placed at 
Rp is Vp{r) = —Vo{a/\r — Rp|), where Vq = e^/ea and e 
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stands for the static dielectric constant of the semicon- 
ductor. 

In order to characterize the electronic eigenstates in the 
impurity band from the point of view of their extended 
or localized nature, we solve numerically the eigenvalue 
problem defined by the Hamiltonian ([T]) for given config- 
urations in which N impurities are randomly placed in a 
three-dimensional volume. Performing the impurity av- 
erage we obtain (Fig. 1) the density of states (DOS) and 
the inverse participation ratio (IPR) for three densities 
on the metallic side of the transition. The impurity band 
develops around the E = level of the isolated impu- 
rity in an asymmetric fashion: the DOS exhibits a long 
low-energy tail while the high-energy part is bounded 
by i? = 1 (in units of Vb). We verify that the width 
of the impurity band increases with the doping density. 
The highest-energy states correspond to electronic wave 
functions localized on small clusters of impurities. This 
clustering is known to happen in realistic systems due to 
the lack of hard-core repulsion between impurities on the 
scale of a 



mm\. 



Before continuing with the analysis of the numerical 
results obtained from the MT model, we discuss some 
technical features of the model and the difficulties that 
we face in trying to improve upon it. Firstly, we notice 
that the chosen basis set is not orthogonal. In princi- 
ple, we can deal with this issue by writing a generalized 
eigenvalue problem which includes the matrix of orbital 
overlaps [l^ E^l- This procedure results in unphysical 
high-energy states (with E :s> I) that necessitate the 
inclusion of hydrogenic states beyond the Is orbital in 
order to be described properly. However, care must be 
taken since enlarging the basis set leads to the problem of 
overcompleteness. Fortunately, for the properties we are 
interested in, the effects arising from non-orthogonality 
are known to be small for moderate doping densities, and 
that is why we will not consider them in the numerical 
work, thus staying within the original MT model. Fi- 
nally, another drawback of the MT model is that the 
high-energy edge of the impurity band overlaps with the 
conductance band, which starts at Vo/2 (the effective Ry- 
dberg) and is not included in the MT description. As seen 
in Fig. [T] the DOS beyond Vb/2 is always very small, and 
therefore we can ignore the effects that the hybridization 
of the bands would yield in a more complete model. 

The IPR measures the degree of extension of the elec- 
tronic wave function. Large values of the IPR denote in 
general extended states while low values are associated 
with localized states. In particular, an homogeneously 
distributed wave function would have an IPK/N — 1, 
while a localized state would exhibit a vanishing IPK/N 
upon increasing values of N. In our study of spin relax- 
ation it is important to know whether the states near the 
Fermi energy are in the region of extended or localized 
states. 

The determination of the mobility edges from the size 
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FIG. 1. Density of states (DOS, thick line and right scale) and 
inverse participation ratio (IPR, left scale) for three different 
densities on the metallic side of the metal-insulator transi- 
tion, obtained through impurity averaging in the Matsubara- 
Toyozawa model. The solid, dashed and dotted curves of 
IPR/iV are for iV = 2744,4096 and 5832, respectively, and 
the vertical lines indicate the Fermi energy. 



scaling of Fig. [T] is not straightforward. This difficulty 
arises from the heavily structured DOS of the MT model. 
At low energy the small values of the DOS translates into 
a poor statistics for feasible sizes. In the high-energy part 
of the impurity band the separation between the curves 
corresponding to different values of N is masked by the 
small values of the IPK/N. 

For the highest density (top panel) the IPR/iV exhibits 
a relatively fiat region at intermediate energies, which 
is approximately independent of N for the two largest 
system sizes. The lower mobility edge can be located 
roughly at E ^ 3.5, where the latter curves separate. 
For lower impurity densities (lower panels) the previous 
analysis becomes increasingly demanding in terms of sys- 
tem sizes. We see that the flat region of IPR/7V shrinks 
from which we can conclude that the lower mobility edge 
is shifting towards higher values of E. 

The study of spin relaxation in doped semiconductors 
with densities close to the that of the MIT calls for a 
generalization of the previously discussed MT model that 
incorporate spin-orbit coupling. Such an extension was 
done in Ref . 15| , where a spin- flip term 



E 



(3) 



was added to T^o (o' 



-a). Similarly to the spin- 
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conserving case, we have 

tZm=Y.^Myp\i^ma). (4) 

The wave function V'mcr is a spin-mixed conduction-band 
state with an envelope part (/>m(r) and a lattice-periodic 
part which has a small spin admixture. In Ref. |15| . the 
expression of the matrix elements of Eq. Q within an 
8-band Kane model were obtained. Since the two-center 
integrals p — m,m' were shown to vanish, the spin-flip 
hopping amplitude is given by the sum of three-center in- 
tegrals with p ^ m,m' . The three-center integrals cannot 
be analytically solved in general. In Ref. [l^ we provided 
approximate analytical expressions of i$^'^,„ which allowed 
us to estimate the spin relaxation times. On the other 
hand, in this work we take the route of the numerical 
evaluation of the three-center integrals. 

Spin-orbit coupling is known to favor the delocaliza- 
tion of disordered systems in two dimensions. In what 
follows we repeat the previous analysis, done for the MT 
model, for the spin-orbit generalized model in this case, 
for increasing values of the spin-orbit coupling Rr (the 
subindex r stands for Rashba). 

For the realistic values of the spin-orbit coupling 
strength {Rr = 1) considered in Ref. [l5|, the spin- 
admixture perturbation energies are, even for largest sys- 
tem sizes that we can treat numerically, orders of mag- 
nitude smaller than the MT level spacing. In Fig. [5] we 
show the spectral decomposition of a MT eigenstate with 
CT = 1 in the basis of spin-admixed eigenstates of Ha+Hi. 
Only for enhanced values of Rr do we obtain significant 
projections into the two spin-admixed subspaces. 

The study of larger values of Rr then appears not only 
as a useful tool for analyzing the progressive inclusion of 
spin-orbit effects, but also as a need for numerical simu- 
lations of the spin dynamics. 

In Fig. El we present the DOS and IPR/iV of the ex- 
tended model for the three densities previously treated 
and various values of the spin-orbit coupling strength Rr . 
The DOS does not noticeably change with Rr, and that 
is why we only present the Rr = 1 case. The spin-orbit 
coupling results in the increase of the IPR/N as a func- 
tion of Rr in the region of extended states. This effect 
is more prominent for the larger density. The low-energy 
sector that has localized states in the MT model exhibits 
IPR/7V curves approximately independent of N, which is 
a signature of the delocalization tendency. 

In Fig. [4] we perform a finite-size scaling of the IPR/iV 
for a given density above the MIT critical density and 
one value of the spin-orbit coupling enhancement factor, 
Rr = 50. 

In conclusion, we revisited the problem of the charac- 
terization of the eigenstates of the Matsubara-Toyozawa 
model from the point of view of their localization, and 
performed a similar analysis in an extended model pro- 
posed recently which includes the structural inversion 
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FIG. 2. Spectral decomposition of a Matsubara-Toyozawa 
eigenstate into the basis set formed by the eigenstates of our 
spin-orbit extended model. The system size is A'^ = 1000 and 
the density is given by (riia^y^^ = 0.33. Top, middle, and 
bottom panels correspond to Rr — 1, 100, 500, respectively. 

asymmetry spin-orbit mechanism. Analyzing the effect 
of spin-orbit coupling of various strengths is necessary in 
order to address the study of the spin dynamics in the 
impurity band of doped semiconductors. We found that 
while the density of states is not considerably modified 
by the spin-orbit interaction, the nature of the states is 
noticeably affected by it showing a tendency to the delo- 
calization. 
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for iV = 2744, 3375 and 4096, respectively, and the vertical 
line indicates the Fermi energy. 
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